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DISTRIBUTED ANTENNA SYSTEM AND METHOD 



ORIGIN OF THE INVENTION 



5 The invention described herein was made by employee(s) of the United States 

government and may be manufactured and used by or for the Government of the United 
States of America for governmental purpose without payment of any royalties thereon or 
therefor. 

1 0 BACKGROUND OF THE INVENTION 



1 . Field of the Invention 

The present invention relates generally to distributed antennas and methods and, 
more specifically, to a distributed antenna with a non-unique phase center that may, in a 
15 preferred embodiment, be used for position and attitude determination. 



2. Background of the Invention 

It is well known that the global positioning system (GPS) may be used for position 
and attitude determinations whereby the GPS receivers solve the pseudo-range equations. It 
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is also well known that much more accurate position and attitude determinations, within less 
than one centimeter, can be made utilizing carrier phase measurements in conjunction with 
the global positioning system (GPS). However, prior art carrier phase location systems may 
frequently experience problems due to their antenna systems. For instance, especially on 
5 moveable items such as vehicles, planes, satellites, thrusters, cell phones, hand-held GPS 
detectors, and the like, the coverage of the antenna(s) may be limited. Rotation or movement 
of the item to which the antenna(s) are secured may result in a loss of signal with respect to 
one or more GPS satellites. While multiple antennas can be utilized to maintain better 
contact, the need for RF switches, lengthy "times to fix", and complex switching algorithms, 

10 may reduce system reliability and performance while significantly increasing the cost 
thereof. Moreover, such antenna systems may not be suitably conformable to the surface of 
the structure upon which they are used so as to be practical, especially for smaller devices 
such as thrusters, cell phones, satellites, and the like. 

While wrap-around antennas, such as spherical and cylindrical antennas, are known 

15 to have a wide angle coverage which would theoretically maintain better contact with a 
plurality of satellites, they also have a non-unique phase center. For instance, if a plurality of 
satellite signals are received by a wrap-around antenna, there may be a plurality of different 
phase centers, equal to the number of observed satellites. The position derived using prior 
art GPS carrier phase techniques requires the position of the antenna phase center. Therefore, 
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the GPS carrier phase solution will be ambiguous when the antenna does not possess a 
unique phase center. The inventor has determined that it would be desirable to remove the 
ambiguity resulting from GPS carrier phase measurements obtained with distributed 
antennas for which the phase center is not unique so as to take advantage of the wide angle 
5 of reception afforded by such antennas. 

Various inventors have attempted to solve problems as indicated by the following 
U.S. Patents. 

U.S. Patent No. 5,890,091 , issued March 30, 1999, to Talbot et al., discloses a fixed 
and roving pair of four observable GPS receivers and a communication link between them 

10 for double differencing code and carrier measurements. Carrier phase integer ambiguities 
are resolved efficiently by searching the simultaneous narrow-lane intersections of both the 
LI and L2 wave fronts propagated by the GPS satellites being tracked. External constraint 
information, such as elevation, is additionally used to speed up integer ambiguity resolution. 
Data between the reference station and the rover is communicated in compressed form at a 

1 5 regular interval, e.g., once a second at each epoch, and demi-measurements of carrier phase 
are obtained more frequently, e.g., ten times a second, and used to propagate solutions 
between epochs. 

U.S. Patent No. 5,943,008, issued August 24, 1999, to D. L. Van Dusseldorp, 
discloses a platform attitude determination made with a single receiver in a global 
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positioning system (GPS). The receiver receives at least three sets of GPS signals from three 
antennas. Two of the GPS signals are delayed so that the receiver receives each GPS signal 
in separate time domain slots. In this way, the synchronization of separate GPS receivers 
does not have to be accomplished. 
5 U.S. Patent No. 5,347,286, issued September 13, 1994, to D. Babitch, discloses a 

system for automatically pointing a directional antenna. The system comprises two GPS 
antennas mounted at horizontally opposed extremities of the directional antenna. The 
placement is such that the GPS antennas lie on a line having a normal vector approximately 
parallel to a boresight of the directional antenna. The outputs from the GPS antennas are 

1 0 down converted to a measurement frequency and differentially phase compared for each of a 
plurality of radio visible GPS satellites. A measurement controller receives the phase 
comparisons. A navigation computer receives measurements from the measurement 
controller and uses a microprocessor to calculate the GPS latitude and longitude and the 
attitude angles of the pair of GPS antennas. A servo points the directional antenna at a 

15 particular target communications satellite, based on the computed azimuth, elevation, 
latitude, and longitude, as provided by the navigation computer. One of these servos is used 
to move the GPS antennas in roughly a circle to average out multipath effects to improve 
accuracy and reduce integer ambiguities. 

U.S. Patent Application No. 5,185, 610, issued February 9, 1993, to Ward et al., 
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discloses a GPS single-receiver pointing/attitude system which derives pointing/attitude 
measurements by correlating a selected GPS code (either P or C/A), recovered from GPS 
navigation signals using a single GPS receiver with multiple GPS antennas ( a reference 
antenna and at least two slave antennas for pointing or three for attitude). For a two antenna 
5 pointing application, with a GPS receiver for each receiver channel, the incoming GPS 
signals are applied to three code correlators assigned to the reference antenna, and three code 
correlators assigned to the slave antenna, which provide corresponding reference and slave I 
and Q correlation outputs. The single-receiver pointing technique involves: (a) using the 
reference I and Q correlation outputs to establish a conventional reference antenna tracking 

10 loop; and (b) processing the reference and slave I and Q correlation outputs (using 
differential carrier doppler phase or code phase measurements) to determine phase 
differences from which pointing can be computed. 

U.S. Patent No. 5,477,458, issued December 19, 1995, to P. V. W. Loomis, discloses 
method and apparatus for providing GPS pseudo range correction information over a 

1 5 selected geographic region S with a diameter of up to 3 OQQ km with an associated inaccuracy 
no greater than 5 cm. N spaced apart GPS fiducial stations, whose location coordinates are 
fixed and known with high accuracy, are provided within or adjacent to the region R. Each 
fiducial station receives GPS signals from at least four common-view GPS satellites, 
compares these coordinates with its known location coordinates, determines the pseudo 
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range corrections for its GPS-determined location, and transmits these correction signals to a 
central station located within or adjacent to the region S. The central station retransmits the 
pseudo range correction signals throughout the region S. A mobile GPS station within or 
adjacent to the region S has stored within it coordinates of the GPS determined last location 
5 of that mobile station and the spatial coordinates of GPS fiducial stations within S that are 
closest to the last determined location of that mobile station. The mobile station then 
computes the differential GPS corrections for the GPS-determined present location of that 
mobile station. Alternatively, the fiducial stations can transmit to the central station 
unprocessed GPS signals for determination of the pseudo range correction signals at the 

10 central station. This approach can be modified if the region S is two-dimensional, where 
only one coordinate is needed. 

U.S. Patent No. 5,587,904, issued December 24, 1 996, to Ben-Yair et al, discloses an 
air combat monitoring system comprising a plurality of GPS receivers, mountable on a 
corresponding plurality of aircraft, an avionics monitor for monitoring the avionics system of 

15 the plurality of aircraft, a memory for storing information received from the GPS receivers 
indicating the location of the plurality of aircraft and information received from the avionics 
monitor in order to provide a reviewable output indication of the performance of the plurality 
of aircraft and an information disseminator for communicating information among said 
plurality of aircraft. 
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U.S. Patent No. 6,018,315, issued January 25, 2000, to Ince et al., discloses a space 
vehicle which uses a signal received from a GPS satellite to determine the attitude of a space 
vehicle. The signal is received through a pair of antennas and converted to sum and 
difference outputs by a sum and difference circuit. These outputs are then conveyed to down 
5 converter/correlators, and then to envelope detectors. The outputs of the envelope detectors 
are then converted to a digital format by an analog to digital converter. The digitally 
formatted sum and difference amplitudes are then used by a processor to calculate the 
required attitude adjustment to the space vehicle. 

U.S. Patent No. 6,052,647, issued April 18, 2000, to Parkinson et al., discloses an 

1 0 automatic control system for land (and possible marine) vehicles based on carrier phase 
differential GPS (CPGPS). The system relies on CPGPS to determine vehicle position 
and attitude very precisely (position to within 1 cm and attitude to within 0.1 degree). A 
system incorporates a technique to calculate and compensate for antenna motion due to 
vehicle roll and pitch. One aspect of the system utilizes an intelligent vehicle controller 

15 that recognizes and adapts to changing conditions, such as vehicle speed, implements 
towed by the vehicle, soil conditions, and disturbance level. The system provides the 
capability to control the vehicle on various paths, including straight lines and arbitrary 
curves. Also described is a technique for initialization and vehicle control using only a 
single pseudolite. 
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U.S. Patent No. 6,061,632, issued May 9, 2000, to S. D. Dreier, discloses a 
system and method for determining which of a plurality of correction signals is to be used 
in a position determination system for determining position. An antenna receives 
broadcasts of correction data on designated frequencies. The received correction data is 
5 analyzed to determine the best correction data signal to use to accurately determine 
position. The correction data from the selected best source of correction data is coupled 
to a position determination system for accurate determination of position. Thus, the 
present invention provides a seamless correction system that determines the best source 
of correction data from available sources and uses the selected correction data to 

1 0 determine location. The seamless correction system updates the source of correction data 
automatically. Thus, the user does not have to take the time and effort necessary to 
determine the location, frequency and format of correction signals broadcast in each area 
and the user does not have to change frequencies or change equipment to obtain accurate 
position determination. 

15 U.S. Patent No. 6,067,484, issued May 23, 2000, to Rowson et al., discloses a 

differential GPS landing system having at least three GPS receivers at known locations in 
spatial proximity to each other, wherein each GPS receiver independently receives GPS 
signals from the plurality of GPS satellites and at least three reference stations, wherein 
each reference station receives a signal from a different two of the at least three GPS 
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receivers and calculates two separate differential corrections for each satellite, wherein 
each differential correction is independently calculated using the signals received from a 
different one of the two GPS receivers. The reference stations preferably average the two 
calculated differential corrections for each satellite to produce an averaged differential 
5 correction for each satellite. The system then validates the averaged differential 
correction for each satellite calculated by each reference station using two integrity 
monitors, wherein each integrity monitor receives the averaged differential corrections 
from a different two of the at least three reference stations, and each integrity monitor 
compares, with respect to each GPS satellite, the averaged differential corrections from 

1 0 the two of the at least three reference stations to produce a validated set of differential 
corrections, and broadcasts the averaged differential correction for each satellite 
calculated by one of the three reference stations. The preferred differential GPS landing 
system further includes a method that utilizes current differential correction error 
statistics, prior differential correction error statistics and a Bayesian interval estimation 

15 bounding process to provide an error bound on the broadcast differential corrections. 

U.S. Patent No. 6,088,653, issued July 1 1 , 2000, to Sheikh et al., discloses A 
method for use in vehicle attitude determination includes generating GPS attitude 
solutions for a vehicle using three or more antennas receiving GPS signals from two or 
more space vehicles. An inertial navigation system is initialized by setting the attitude of 
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the inertial navigation system to a GPS attitude solution generated for the vehicle and/or 
the attitude of the inertial navigation system is updated using the GPS attitude solutions 
generated for the vehicle or GPS estimated attitude error generated for the vehicle. A 
system for use in vehicle navigation is also provided. The system generally includes three 
5 or more GPS antenna/receiver sets associated with a vehicle, an inertial measurement unit 
that provides inertial measurements for the vehicle, a processing unit of the system 
having the capability for generating GPS attitude computations for the vehicle using the 
three or more GPS antenna/receiver sets and signals from two or more space vehicles; the 
GPS attitude computations include at least one of absolute attitudes and estimated attitude 

10 errors. The processing unit of the system also includes a filter for generating estimates of 
attitude for the vehicle using the inertial measurements from the inertial measurement 
unit and the attitude computations. 

U.S. Patent No. 6,101,430, issue August 8, 2000, to Fuller et al., discloses an 
improved GPS Attitude Receiver for spacecraft uses concurrent line-bias estimates 

15 derived from GPS signals for attitude determination. Attitude determination is obtained 
from GPS measurements by an interactive solution of the known "Attitude Master 
Equation", using RF differential phase measurements and calibrated parameters of 
baseline vectors and line-biases. A "weighted fit error" W, is derived from the sum 
squared of discrepancies between predicted and actual measurements, weighted by 
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measurement variances. By minimizing W, the line-bias estimations obtained give more 
current data including thermal effects and eliminate separate ground calibration tests. 

The above prior art does not disclose a system that takes advantage of a greatly 
improved wrap around coverage from one or more distributed antennas with non-unique 
5 phase centers. Therefore, those skilled in the art have long sought and will appreciate the 
present invention that addresses these and other problems. 
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SUMMARY OF THE INVENTION 

An object of the present invention is to provide an improved distributed antenna and 
method. 

5 Yet another object of the present invention is to provide an antenna system, that may 

be used with GPS applications or other multi-transmitter systems, which has a wide angle, or 
even completely spherical, field of view without the introduction of error due to the non- 
unique phase center such that as a body supporting the antenna system rotates or partially 
blocks the antenna system field of view with respect to one or more transmitters, contact 
1 0 with all transmitters is continuously maintained. 

Another object of the present invention is to provide an antenna system that is 
conformable to a surface of a body including smaller bodies such as cell phones, small 
satellites, thrusters, and the like. 

One preferred feature of the invention permits continuous operation of a single 
1 5 distributed antenna without switching between antennas. 

Another feature of a presently preferred embodiment is a conformable antenna for 
moveable spherical or cylindrical bodies. 

An advantage of the present invention is the possibility of elimination of lengthy time 
to first fix following antenna switching. 
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Yet another advantage of the present invention is the ability to reduce or eliminate the 
integer ambiguity problem associated with GPS attitude determination. 

Yet another advantage of the present invention is the possibility of complete attitude 
solution (three degrees of freedom) with only two antennas, or a partial attitude solution (two 
5 degrees of freedom) with a single antenna. 

These and other objects, features, and advantages of the present invention will 
become apparent from the drawings, the descriptions given herein, and the appended claims. 
It will be understood that above listed objects, features, and advantages of the invention are 
intended only as an aid in understanding aspects of the invention, are not intended to limit 
1 0 the invention in any way, and do not form a comprehensive list of such objects, features, and 
advantages. 

Therefore, the present invention comprises a method for utilizing a plurality of 
transmitters to determine one or more location characteristics of a body wherein the plurality 
of transmitters may produce a plurality of RF carrier signals. The method may comprise one 
15 or more steps such as mounting one or more distributed antennas to the body. Each of the 
one or more distributed antennas may have a non-unique phase center. Other steps may 
comprise receiving the plurality of RF carrier signals from the plurality of transmitters with 
the one or more distributed antennas, and determining the one or more location 
characteristics of the body. In a preferred embodiment, the method further comprises 
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utilizing carrier phase measurements for determining the one or more location characteristics 
of the body. The one or more location characteristics may comprise a position and/or an 
attitude of the body. In one embodiment of the invention, each of the one or more 
distributed antennas has a substantially spherical coverage. The invention may be 
5 successfully utilized even when at least a portion of a view of the one or more transmitters 
by the one or distributed antennas is blocked by the body. In one embodiment, each of the 
one or more distributed antennas are circularly constructed to thereby have a respective 
physical origin. The method may further comprise determining values related to a vector 
from the respective physical origin of the body or antenna to a phase center with respect to 

1 0 each of the plurality of transmitters. 

Other steps may comprise determining a distance between each of the plurality of 
transmitters and each of a plurality of phase centers whereby each phase center corresponds 
to one of the plurality of transmitters. In one embodiment, the method comprises obtaining a 
coarse estimate of a position of the body, and utilizing the coarse estimate for determining a 

15 unit vector related to a first position vector of a physical center of the distributed antenna 
with respect to a reference system. Continuing, the method may further comprise utilizing 
the unit vector so obtained for determining an improved estimate of a position of the body, 
and utilizing the improved estimate for determining the unit vector more accurately. 
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The method may also comprise determining information related to a phase center 
constellation comprised of a plurality of phase centers such that each phase center in the 
phase center constellation is related to a respective of the plurality of transmitters. In one 
embodiment, the method comprises determining a plurality of vector magnitudes whereby 
5 each vector magnitude is related to a vector from a physical center of the distributed antenna 
to one of the plurality of phase centers. 

In other words, the method may comprise utilizing a plurality of transmitters in a 
plurality of locations to determine one or more location characteristics of a body including 
steps such as mounting one or more distributed antennas to the body wherein each of the one 
10 or more distributed antennas may have a plurality of phase centers with respect to the 
plurality of locations of the plurality of transmitters, and determining a plurality of values 
related to the plurality of phase centers. 

Specific embodiments may include determining a vector to a reference center of the 
body from a fixed coordinate system and/or determining a vector from the reference center to 
15 an i th satellite and/or determining a vector from the reference center to a phase center related 
to the i th satellite. Other steps may include determining an attitude vector for the body and/or 
determining a unit vector for the body from the reference center with respect to an i th 
satellite. In one embodiment, a step involves estimating the unit vector by determining a 
position vector of the body with respect to a coordinate system. If desired, a further step 
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may comprise reducing the error of the estimate of the unit vector by iteration techniques. 

In one embodiment, the method comprises steps such as measuring a carrier phase 
from an i th satellite, and adding a correction to the measured carrier phase. The correction 
may be determined utilizing a known attitude of the body. In another embodiment, steps 
5 may include determining a magnitude of a vector from the reference center to a phase center 
related to the i th satellite. 

The invention may comprise a system for determining location characteristics of a 
body utilizing a plurality of spaced apart transmitters. In one embodiment, the system may 
comprise one or more antennas mounted to the body wherein each of the one or more 
10 antennas may have a non-unique phase center with respect to the plurality of spaced apart 
transmitters. The one or more antennas may preferably provide a wide angle coverage for 
maintaining contact with the one or more spaced apart transmitters. As well, means are 
provided for determining the location characteristics in response to reception of signals from 
the spaced apart transmitters by the one or more antennas. The means for determining the 
1 5 location characteristics may comprise carrier phase related equations. 

The invention may comprise a radiator system for determining location 
characteristics of a body utilizing a plurality of spaced apart transmitters wherein the body 
may have a curved surface. The radiator system may then comprise one or more radiators 
mounted to the curved surface of the body so as to conform to the curved surface. Each of 
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the one or more radiators may have a non-unique phase center with respect to the plurality of 
spaced apart transmitters and the one or more antennas providing a wide angle coverage for 
maintaining contact with the one or more spaced apart transmitters. In one embodiment, the 
one or more radiators comprise a circular ring. In another embodiment, the one or more 
5 radiators comprise a plurality of circular rings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



FIG. 1 is a schematic showing a wrap-around antenna on a spherical body in accord 
with the present invention, an i th satellite, and an associated i th antenna phase center; 
5 FIG. 2 is a schematic showing an orientation of a body z-axis in fixed coordinates in 

accord with the present invention; 

FIG. 3 is a schematic showing local coordinates of the spherical body of FIG. 1 in 
accord with the present invention; 

FIG. 4 is a graph showing antenna phase variation over an elevation cut and sine 
10 function; 

FIG. 5 is a schematic showing two wrap around antennas mounted on a cylindrical 
structure in accord with the present invention; and 

FIG. 6 is a schematic view of a spherical antenna pattern of a distributed antenna in 
accord with the present invention which may be utilized for maintaining signal contact with 
15 multiple transmitters as a body rotates without the need for RF switches. 

While the present invention will be described in connection with presently preferred 
embodiments, it will be understood that it is not intended to limit the invention to those 
embodiments. On the contrary, it is intended to cover all alternatives, modifications, and 
equivalents included within the spirit of the invention and as defined in the appended claims. 
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DESCRIPTION OF PREFERRED EMBODIMENTS 



The present invention provides apparatus and methods for a distributed antenna 
utilized for the determination of the position and/or attitude of a body when the observable 
5 is the phase of the RF carrier signal. In one presently preferred mode of operation, the 
distributed antenna receives signals from the constellation of the global positioning system 
(GPS) satellites. However, other transmitters could be used either in conjunction with GPS 
satellites or instead of GPS satellites. The phrase "distributed antenna" as used herein refers 
to an electromagnetic wave radiator with a non-unique phase center, i.e., the apparent phase 
1 0 center is dependent upon the relative position of the observer (or GPS satellite). A radiator 
in accord with the present invention may comprise a single, continuous radiating element or 
may comprise two or more discrete radiating elements whose received signals are 
combined. 

The antenna of the present invention is preferably designed to provide maximum gain 
15 over a desired solid angle in which coverage is required. In one preferred embodiment, 
nearly complete spherical coverage, i.e., 4n steradians, is provided to the extent possible 
without introducing RF switches for antenna selection. For one example of a suitable 
radiation pattern, FIG. 6 discloses an antenna with spherical radiation pattern 1 1 in accord 
with the present invention. It is commonly known that the pattern of the "wrap-around" 
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antenna has nulls (not shown) in the direction of the + and - z-axis (i.e., theta = 0 and theta 
= 1 80 degrees.) The width of the nulls tends to zero degrees as the body diameter increases. 
Thus, a rotatable body located at the origin of radiation pattern 1 1 may remain in contact 
with multiple satellites even as the body is rotated without the need for RF switches. 
5 Antennas commonly used in GPS determination of position and attitude have a 

unique phase center that is known a priori. The GPS receivers may solve the well-known 
pseudo range equations (Eqn. 1) for the three unknowns (x, y, and z) that describe the 
position of the body, as well as the time bias of the user's GPS receiver. 



1 0 \\r, -r B \\ = v,. + ct h , / = 1,2,3,4 Eqn. 1 

where: 

r si : position of the /"' satellite 

r n : position of the receiver 

v i : pseudo - rangemeasurement 

c : velocity of light 

t h : user clock bias 

The accuracy of the position solution may be enhanced by using GPS carrier phase 
1 5 measurements rather than the pseudo-range measurements. However, the position derived 
using the GPS carrier phase utilizes the position of the antenna phase center. Therefore, the 
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GPS carrier phase solution will be ambiguous when the antenna does not possess a unique 
phase center. 

One intention of the method described herein is to remove the ambiguity resulting 
from GPS carrier phase measurements obtained with a distributed antenna for which the 
5 phase center is not unique. 

Referring now to the drawings, and more particularly to FIG. 1, body 10, whose 
position/attitude is to be monitored, may for instance be a metallic sphere such as might be 
used for a satellite. However, other shapes for body 1 0 could also be used such as cylinders 
and/or a suitable shaped antenna conformable to the surface of the body whose 

10 position/attitude is to be monitored. Again for purposes of explanation, antenna 12 will be 
assumed a wrap-around or circumferential belt type of antenna as depicted in FIG. 1 . It will 
be understood that wrap-around antenna 1 2 is intended to be representative of a wide variety 
of antennas that can be placed in a circular configuration or in a piecewise linear 
approximation to a circular configuration. The method described herein is applicable even 

15 when the body internal to the wrap-around antenna, such as body 10, obscures the view or 
visibility of part of the antenna, i.e., some part of the antenna or some antennas are not 
significantly contributing to the composite signal due to blockage of the line of sight by the 
body. While body 1 0 is a sphere, the same reasoning would also apply to a cylinder and/or 
any distributed antenna or distributed antenna system possessing suitable symmetry. 
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In FIG. 1 5 position vector 14 to center 16, within fixed coordinate reference system 
18, is shown as r b . The z-axis 20 of sphere 10, which is normal to the plane containing 
wrap-around antenna 12, is denoted by z B . The position of the i th satellite, noted as satellite 
24, with respect to the origin of reference system 18 is vector r si (not shown). Vector 26 is 
5 then r si - r b . Vector 22, labeled as vector r p j, is the vector from center 1 6 of sphere 1 0 to the 
phase center as seen by the i th satellite, e.g., satellite 24. Due to the symmetry of wrap- 
around antenna 12, vector 22, i.e., vector r P j, lies at the intersection of the plane formed by 
vectors 20, 26 and the equatorial plane containing wrap-around antenna 1 2. Although vector 
22, i.e., vector r p j, is shown in FIG. 1 to have a magnitude equal to the radius of sphere 10, 
1 0 this is not generally the case. However, as shown subsequently, the magnitude of vector 22, 
i.e., vector r P j, can be estimated, and is essentially constant as long as the i th satellite, such as 
satellite 24, is in the same plane. 

When GPS carrier phase measurements are employed, the magnitude of the vector 28 
(or the " || Vj || = v" in Eqn. 1 ) is the observable and can be considered to be the distance from 
15 the phase center of wrap-around antenna 12 to the i th satellite, such as satellite 24, as 
determined by the carrier phase measurement. 

In a presently preferred embodiment of the present invention, the following equations 
(Eqn. 2) are hereby proposed to describe the case when antenna 12 is a wrap-around or 
distributed antenna. 
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IHf = \K ~r B f +\\r pi f -2r pi .(r s , -r B ) 
or, equivalently, 

5 

l v /| 2 =h -r B ( +\\r pi f -2\\r pi l\r xi -r B \\cos(/},) Eqn. 2 

* 

where 

cos(/?,) = sin(a,) = ^ 

10 and 

z B = sm(9 B \x, cosfo, )+ y e sin(^ )]+ z e cos{0„ ) 

Eqn. 2 is derived from the "law of cosines" and trigonometric relationships. 9 B and 
15 <f> B describe the attitude of vector 30 shown in FIG. 2, designated as z B , for sphere 1 0 of FIG. 
1 and FIG. 3. Eqn. 2 is a non-linear set of equations with five unknowns (assuming 0 B and 
<|)b are not already known) or three unknowns (if 9 B and <|>b are already known) and a number 
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of equations equal to the number of observed satellites, such as satellite 24. If the user 
desires position, the primary unknowns are the three vector components of r b . If the clock 
bias is also corrected, an additional unknown is necessary. It should be noted that the 
attitude of the body relative to the fixed coordinates might be described in coordinates other 
5 than those shown in FIG. 2, i.e., 0b and fa. 

FIG. 3 shows the local coordinates of sphere 10. Note that the rotation of sphere 10 
about the axis zb, i.e., the location of local body axes xb and ye, is not required for the 
solution, nor can it be determined by this method when the antenna has complete azimuthal 
symmetry. Furthermore, the attitude of vector 30, or zb, is ambiguous in sign when Eqn. 2 is 
1 0 used to determine 0b and <|>b. This, however, does not impact the position solution. 

As an alternative to the set of equations of Eqn. 2, it may be recognized that the term 

^r pi ^ is much smaller than the remaining 1 st or 2 nd equations in Eqn. 2. Recognition of this 
can be used to simplify the 1 st equation in Eqn. 2 to the following Eqn. 3. 



15 




or, equivalently, Eqn.3 
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IklhlK-^l-Hcos^,.) 

where Pi is defined, as in Eqn. 2, : 



cos (A ) = sin («,)= 1- 



(f si -f B ) 



r si ~r B 



Eqn. 4 



where r sjh = ^ — ^) 



Ki ~ r B \\ 



1 0 Note that r sib is a unit vector. Furthermore, since \r sj || » \r h || (normally true for GPS), the 
following approximation of the position vector can be made for the purposes of obtaining the 
unit vector r : 



15 



Eqn. 5 
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where r h is a coarse estimate of the body or sphere 10 position. This coarse estimate may 
come from various sources such as, for instance: 
5 i) A position solution using wrap-around antenna 1 2 with no phase corrections; 

ii) A position solution using wrap-around antenna 12 with the pseudo-range 
equations; or 

iii) A previous solution of the carrier-phase range equations with corrections. 

The error in r sib should be negligible for most GPS application. However, if desired, the 
10 error may be further reduced by iterating. That is, Eqn. 2 or Eqn. 3 can be solved the first 
time using the approximation of Eqn 5. Then, the first solution, 7 h] , can be used in Eqn. 5 

for a revised estimate of r sih which can them be used to obtain a second solution, 7 h2 , and so 
forth, to the degree of accuracy necessary. 

In another embodiment, the method of the present invention may be applied to the 
15 carrier phase range equations to solve for the position of the body when the attitude of the 
body, such as sphere 10, is known a priori. The adaptation allows the carrier phase range 
equations to be solved using existing algorithms. The carrier phase measured from each 
satellite is modified by a correction. The correction is the final term in Eqn. 3. This term 
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may be found using the approximation r h for rb as describe above. Since the attitude of the 
body is known, the body z-axis is known, and the final term in Eqn. 3 can be computed. 
This term is then added to the measured carrier phase from the i th satellite: 



The resultant, then represents the carrier-phase from the i satellite that would be measured 
at the origin of the body. 

The present invention permits determination of the phase center constellation for 

10 wrap-around antenna 12. As described above, the vector r p j can be assumed to lie at the 
intersection of two planes. One plane is the plane spanned by the body z-axis 20, zb, and the 
vector to the i th satellite, r S i. The second plane is the plane containing the wrap-around 
antenna, such as wrap-around antenna 12. This assumption is base on the symmetry of the 
radiating structure. It remains only to determine the magnitude of vector 22, r p j. If the phase 

15 is measured over a plane cut with constant azimuth angle in the local body coordinates, i.e., 
\|/ is constant in FIG. 3 and £ varies, then the phase at the phase center will lead the phase at 
the local coordinate origin 16, approximately, by: 



K/- r *ll = ll v /ll + ll r /rf|l cos to) 



Eqn. 6 



A = *|rjsin(£) 



Eqn. 7 
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Therefore, by measuring the phase of the wrap-around antenna versus with v|/ constant, the 
magnitude of the vector r P i can be determined by a standard least-squares or other appropriate 
fitting technique applied to Eqn. 8: 

5 

Aj =a 0 + k\\r pi \\sm(Z),j = l,...,N Eqn. 8 



where oto is an additive constant, and the index j ranges over the N sample points. A 
multitude of cuts, each at different angles, \\t 9 and each corresponding to a different value for 
1 0 the index i, can be used to determine an average value for the magnitude of the vector 22, r pi . 
FIG. 4 shows the simulated phase variation as the angle £ is varied between 0° and 90° as 

well as a function given by the right hand side of Eqn. 3 with chose to satisfy Eqn. 3 

exactly at £= 90°. 

Therefore, the foregoing disclosure and description of the invention are illustrative 
15 and explanatory thereof, and various changes in the method steps and also the details of the 
apparatus may be made within the scope of the appended claims without departing from the 
spirit of the invention. 
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